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I.   INTRODUCTION 

Wave  theories  developed  for  deep  water  waves  can  be 
applied  with  some  degree  of  certainty,  and  can  be  tested 
in  laboratory  and  field  situations.   These  theories  general 
ly  hold  until  such  time  that  the  shoaling  process  begins, 
and  with  somewhat  lesser  accuracy,  throughout  the  shoaling 
process  up  to  the  point  of  near  breaking.   At  the  breaker 
point,  however,  there  is  a  transition  from  ordered  to 
turbulent  motion  and  the  description  of  wave  kinematics 
becomes  more  difficult.   The  most  forthright  approach  to 
the  problem  of  describing  the  kinematics  of  surf  zone 
breaking  waves  is  through  direct  measurements.   Advances  in 
instrument  design  have  led  to  simple,  sturdy,  measuring 
devices  with  rapid  response  time  which  can  sense  the  small 
scale  as  well  as  the  large  scale  features. 

The'  study  of  the  kinematics  of  breaking  waves  in  the 
surf  zone  has  progressed  slowly  due  both  to  the  problems 
encountered  in  making  direct  field  measurements  and  the 
difficulty  in  modeling  the  surf  zone  in  the  laboratory. 
Breaking  waves  were  measured  in  the  laboratory  in  early 
work  by  Iverson  (1953)  using  photographic  techniques. 
The  length  of  the  channel  limited  the  wave  type  to  plunging 
and  surging  breakers.   Adeyemo  (1970)  made  similar  labora- 
tory measurements  using  hydrogen  bubble  and  photographic 
methods.   Gaughan  and  Komar  (1975)  applied  the  theory  of 
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wave  propagation  in  water  of  gradually  varying  depth,  as 
developed  by  Biesel,  to  determine  the  dependence  of  breaker 
type  on  the  beach  slope  tangent  and  the  deep  water  wave 
steepness . 

Inman  and  Nasu  (1956)  made  field  measurements  of  water 
particle  velocity  by  measuring  the  drag  force  under  the 
wave  in  order  to  infer  particle  velocities.   Miller  and 
Ziegler  (1964)  used  both  acoustic  and  electromagnetic 
current  meters  to  determine  the  particle  motion  in  the 
surf  zone.   Walker  (1969)  made  studies  using  propeller- type 
flow  meters.   Wood  (1973)  measured  waves  and  currents  in 
the  surf  zone  using  movies  of  dye  movement  and  capacitance 
wave  gauges.   Fiihrboter  and  Biisching  (1974)  utilized  a  two- 
component  current  meter  and  two  pressure- type  wave  meters 
to  measure  simultaneous  orbital  velocities  and  water  levels 
Thornton  (1968),  Steer  (1972),  Thornton  and  Richardson 
(1973)  and  Bub  (1974)  used  pressure  meters,  capacitance 
wave  gauges  and  electromagnetic  current  meters  to  measure 
surface  profiles  and  particle  velocities;  the  work 
presented  here  is  an  extension  of  these  studies. 

The  objective  of  this  research  is  to  study  the  kine- 
matics of  water  particles  in  breaking  waves  within  the 
surf  zone.   Simultaneous  measurements  were  made  of  the 
instantaneous  wave  profile  and  the  horizontal  water 
particle  velocities  at  fixed  locations  in  the  surf  zone. 
Estimates  of  the  probability  density  functions  and  power 
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spectra  of  the  wave  heights  and  particle  velocities  are 
made.   The  applicability  of  using  linear  wave  theory  as 
a  spectral  transfer  function  in  computing  velocity  spectral 
components  from  the  power  spectrum  of  the  waves  is 
measured.   The  computed  velocity  spectra  are  compared  with 
actually  measured  velocity  spectra.   Theoretical  phase 
spectra  between  wave  gauges  are  calculated  and  compared 
to  measured  phase  spectra. 
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II.   MEASUREMENTS 

A.   EXPERIMENT  SITES 

The  experimental  sites  are  in  the  vicinity  of  Monterey, 
California.   The  beaches  here  are  some  of  the  first  inten- 
sively studied  to  gain  an  understanding  of  amphibious  war- 
fare techniques  and  are  described  by  Bascom  (1964) .   It  was 
desired  to  measure  the  various  types  of  breaking  waves  in 
the  study  including  spilling,  plunging  and  surging  breakers. 
The  manner  in  which  waves  break  depends  very  much  on  the 
characteristics  of  the  beach  and  near-shore  bottom  slope 
(Table  I) .   Plunging  and  spilling  breakers  occur  most 
frequently  at  the  Del  Monte  Beach  site  within  Monterey  Bay. 
The  waves  at  this  location  are  generally  severely  direc- 
tionally  filtered  and  refracted  by  the  geometry  of  the  bay 
and  impinge  almost  perpendicular  to  the  shore  with  a 
resulting  simplification  in  the  wave  description.   A  second 
experimental  site  was  Carmel  River  Beach,  five  miles  to 
the  south,  where  the  beach  is  very  steep  and  very  often 
has  surging  type  breakers.   Again  here,  the  beach  is  within 
an  embayment  and  the  waves  impinge  almost  perpendicular  to 
the  shore.   A  typical  beach  profile  and  instrument  loca- 
tion for  Del  Monte  Beach  is  shown  in  Figure  1.   Del  Monte 
Beach  profiles  for  4-10  March  1975  are  given  in 
Appendix  A.   Figure  2  shows  the  beach  profile  and  instrument 
location  at  Carmel  River  Beach  on  29  May  1975.   A  beach 
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profile  taken  the  following  day  was  identical  to 
Figure  2. 

B.   INSTRUMENTATION 

Many  problems  were  encountered  in  the  initial  attempt 
in  February  1975  to  measure  the  desired  parameters.   The 
first  wave  gauges  were  manufactured  from  1.0  cm  diameter 
stripped  RG-11  coaxial  cable  as  was  done  by  Bub  (1974) . 
Large  masses  of  kelp  torn  loose  by  high  waves,  flotsom 
and  the  pounding  surf  easily  broke  these  cables.   Scouring 
around  the  base-plates  of  the  instrument  towers  and  the 
blade-type  anchors  caused  the  towers  to  topple.   The  con- 
struction of  the  wave  gauges  was  modified  to  a  three- 
eighths  inch  outside  diameter  stainless  steel  rod,  covered 
with  a  one-sixteenth  inch  wall  thickness  teflon  tubing 
fitted  tightly  over  the  rod;  these  wave  gauges  withstood 
the  severe  forces  in  the  surf  well.   Steel  pipe  extensions 
each  0.6  m  in  length  were  added  under  the  tower  base  plates 
and  resulted  in  stable  platforms.   An  extension  was  built 
for  the  top  of  the  anchors;  with  this  device,  the  anchors 
were  sledge-hammered  deeper  into  the  beach  so  that  the  top 
of  each  anchor  was  approximately  15  cm  below  the  sand.   The 
extension  was  then  removed  and  scouring  was  minimized. 

Measurements  at  Del  Monte  Beach  were  made  using  two 
electromagnetic  flow  meters  and  three  capacitance  wave 
gauges.   One  flow  meter  and  two  wave  gauges  were  used  at 
Carmel  River  Beach.   The  instruments  were  mounted  on  the 
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towers  within  the  surf  zone.   All  equipment  was  calibrated 
in  the  laboratory  prior  to  the  experiments. 

The  electronics  package  for  each  of  the  2.5  m 
capacitance  wave  gauges  was  constructed  using  the  design 
of  McGoldrick  (1969) .   The  capacitance  wave  gauge  operates 
on  the  principle  that  a  change  of  capacitor  plate  dimensions 
changes  the  circuitry  voltage  output.   The  rod  and  the 
sea  water  act  as  plates  and  the  teflon  insulation  as  the 
dielectric.   As  the  wave  elevation  fluctuates,  the  capaci- 
tance of  the  circuit  changes  and  the  voltage  output  responds 
linearly.   These  water  surface  fluctuations  were  sensed  by 
a  transistorized  circuit  operated  by  self-contained  bat- 
teries.  The  electronics  package  was  housed  in  a  water- 
tight brass  case  which  was  mounted  on  the  tower.   This 
enabled  the  connecting  leads  to  be  less  than  30  cm,  thereby 
minimizing  wire-to-wire  capacitance.   Accuracy  was  estimated 
to  be  ±0.005  m.   The  calibration  plots  are  shown  in 
Appendix  B.   Minor  modifications  to  the  wave  gauges 
immediately  following  the  Del  Monte  Beach  measurements 
are  the  cause  for  the  different  calibration  graphs  in 
Appendix  B. 

The  two  flow  meters  were  Marsh-McBirney  Model  721 
Electromagnetic  Current  Meters.   The  flow  meter  operation 
is  based  on  Faraday's  principle  of  electromagnetic  induction 
Each  probe  measures  water  velocity  in  two  orthogonal 
directions.   The  flow  meters  were  calibrated  with  an 


oscillating  platform  attached  to  an  eccentric  arm  driven 
by  a  variable  speed  motor  using  the  method  of  Thornton 
and  Krapohl  (1974) .   Measurement  accuracy  was  determined 
to  be  ±0.005  m/sec  during  calibration.   The  flow  meter 
calibration  plots  are  shown  in  Appendix  C. 

In  the  field  the  instruments  were  attached  to  6.3  cm 
outside  diameter  steel  pipes  which  were  3.6m  high  with  a 
1.0  m  base  plate  and  the  0.6  m  steel  pipe  extensions.   The 
towers  were  placed  on  a  line  perpendicular  to  the  predomi- 
nant wave  direction  and  were  erected  during  low  tide  when 
the  beach  was  easily  accessible.   The  measurements  were 
then  conducted  at  high  tide.   The  tidal  range  in  the 
Monterey  area  is  typically  two  meters.   The  towers  were 
supported  by  steel  guy  wires  fastened  to  the  anchors. 
Several  types  of  anchors  were  tried  and  the  blade  anchor 
worked  quite  satisfactorily  in  the  relatively  coarse  sand. 
The  wave  gauges  were  attached  to  the  towers  with  locally 
fabricated  micarta  clamps.   A  1.6  cm  steel  stiffening  rod 
was  added  to  the  arrangement  in  order  to  reduce  wave  gauge 
vibration  following  breaker  impact.   The  flow  meters  were 
positioned  directly  under  the  wave  gauges  with  the  axes 
aligned  horizontally  (+X,  shoreward)  and  vertically  (+Y, 
upward).   A  carpenter's  level  was  used  to  establish  axis 
alignment  with  an  estimated  error  of  ±2  degrees.   A  typical 
tower  and  sensor  arrangement  is  shown  in  Figure  3. 

All  signals  were  cabled  ashore  and  recorded  on  a  Vidar 
Corporation  32-channel  digital  data  acquisition  system. 
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A  Sangamo  model  3500  14-channel  FM  tape  recorder  was 
utilized  as  a  secondary  recording  system.   A  Brush  8-channel 
strip  chart  recorder  was  used  to  monitor  the  instrumentation 
performance  during  recording  and  as  a  means  to  select  the 
appropriate  data  sections  to  be  analyzed. 
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III.   ANALYSIS  OF  DATA 

Record  lengths  of  30  minutes  from  each  data  set  were 
analyzed.   The  details  of  the  analysis  system  including 
calibration  factors  are  given  in  Appendices  D  and  E.   A 
mean  value  was  calculated  for  all  data  sets  and  the  data 
was  linearly  detrended  to  exclude  the  rise  and  fall  of  the 
tide.   The  variance,  standard  deviation  and  average  period 
were  calculated.   The  average  period  was  determined  by 
calculating  the  time  between  zero  upcrossings.   However, 
perturbations  caused  by  secondary  gravity  waves  "riding" 
the  primary  gravity  waves,  capillary  waves  and  minor  instru- 
ment noises  increased  zero  upcrossing  occurrences  and 
caused  the  calculated  average  periods  to  be  lower  than 
visually  observed.   Probability  density  functions  for  each 
data  set  were  calculated  and  graphically  compared  with 
Gaussian  and  Gram-Charlier  distributions. 

For  each  signal  an  auto-covariance  function  was  calcu- 
lated and  smoothed  with  a  Parzen  window.   A  Fourier  trans- 
form was  then  applied  to  the  smoothed  auto-covariance 
function  and  the  power  spectrum  determined.   A  cross- 
covariance  function  between  data  sets  was  computed,  smoothed 
with  a  Parzen  window,  a  Fourier  transform  applied,  and  the 
cross-spectrum  computed.   The  coherence  and  phase  were  then 
calculated.   To  compare  waves  and  flow  velocity,  the  wave 
profile  spectrum  was  converted  to  a  theoretical  velocity 
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spectrum  for  comparison  with  the  measured  flow  velocity 
spectrum.   The  two  power  spectra,  the  coherence  and  the 
phase  were  then  plotted.   In  comparing  two  wave  profiles, 
theoretical  celerities  were  calculated  using  linear  theory 
and  its  shallow  water  approximation  which  were  then  con- 
verted to  phases  and  plotted  for  comparison  with  the 
measured  phase  difference. 

Aliasing  testing  was  performed  on  all  spectra. 
Nyquist  frequencies  of  1.25  and  0.98  Hz  were  determined 
optimal  for  Del  Monte  Beach  and  Carmel  River  Beach  data 
sets,  respectively.   The  frequency  bands  of  0.0  to  1.25 
and  0.98  Hz  include  the  region  of  gravity  waves  (0.033 
to  1.0  Hz) . 

The  maximum  lag  time  in  calculating  the  covariance 
functions  was  taken  as  five  percent  of  the  record  giving 
a  spectral  bandwidth  resolution  of  0.00556  Hz  and  results 
in  40  degrees  of  freedom  for  each  spectral  estimate.   The 
80  percent  confidence  limits  for  40  degrees  of  freedom 
using  a  chi-square  distribution  are  found  to  be  between 
0.73  and  1.30  of  the  measured  power  spectral  estimates. 
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IV.   RESULTS 

A.  QUALITATIVE  DESCRIPTION 

A  number  of  similarities  of  wave  form  can  be  observed 
for  various  types  of  breakers  occurring  on  different 
beaches  around  the  world.   Figures  4  and  5  are  typical 
analog  records  of  waves  and  horizontal  velocities  beneath 
the  waves  obtained  from  Del  Monte  and  Carmel  River  Beaches, 
respectively.   In  general,  there  is  a  quick  drawdown  of 
water  just  before  the  breaking  wave  arrives,  followed  by 
a  steep,  almost  vertical  leading  edge,  and  a  sloping  pro- 
file toward  the  trailing  edge.   On  the  trailing  edge  of 
the  wave  secondary  waves  are  often  noted.   These  are 
harmonics  of  the  primary  wave  frequency  and  are  indicative 
of  very  non- linear  waves.   At  Del  Monte  Beach  the  waves 
break  as  plunging  and  spilling  breakers  giving  a  generally 
saw-toothed  shape;  the  plunging  or  spilling  occurs  rapidly 
at  the  crest  and  moves  down  the  wave.   The  surging  breakers 
occurring  at  Carmel  River  Beach  rise  and  fall  off  more 
gradually. 

B.  PROBABILITY  DENSITY  FUNCTIONS 

1.   Gaussian  and  Gram-Charlier  Frequency  Distributions 
Frequency  distributions  were  computed  for  each 
measured  quantity  for  all  runs  and  compared  with  Gaussian 
and  Gram-Charlier  distributions.   The  chi-square 
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goodness-of- f it  test  was  computed  using  both  theoretical 
distributions.   The  variance,  standard  deviation,  skewness 
and  kurtosis  were  also  computed  for  all  data  sets  and  are 
summarized  in  Table  II.   The  closer  the  chi-square  fit 
parameter  is  to  zero,  the  better  the  fit. 

In  deep  water  ,the  sea  surface  has  approximately  a 
Gaussian  distribution  which  results  in  a  skewness  of  zero 
and  a  kurtosis  of  3.   However,  non-linearities  introduce 
skewness  and  kurtosis  values  that  deviate  from  Gaussian 
and  result  in  a  distribution  more  closely  approximated  by 
the  Gram-Charlier .   These  non-linearities  are  evident  in 
the  frequency  distribution  of  Wave  Gauge  #2  on  4  March  as 
shown  in  Figure  6,  and  result  in  bi-modal  or  tri-modal 
distributions.   Figure  7  is  the  graph  of  distributions 
computed  for  the  particle  velocity  at  Flow  Meter  #2  on  the 
same  date.   The  latter  has  a  uni-modal  distribution. 
Appendix  F  contains  distributions  computed  for  all  other 
data  sets. 

The  distributions  for  all  flow  velocities  are  either 
uni-modal  or  have  less  pronounced  secondary  peaks  compared 
to  the  distributions  of  wave  heights.   As  noted  in  Table  II, 
this  results  in  nearly  equal  fit  parameters  for  both  Gaussian 
and  Gram-Charlier  distributions,  although  the  latter  values 
are  smaller,  and  therefore  give  the  better  fit  in  ten  of 
eleven  cases.   The  exception  is  Flow  Meter  #2  on  8  March. 
Additionally,  the  fit  parameters  for  velocity  components 
are  much  closer  to  zero  than  are  the  values  computed  for 


27 


fn 

a> 

H 

CD 

(I) 

rH 

rH 

rC 

-P 

fn 

CO 

CD 

CN 

LO 

rH 

r-- 

CO 

LO 

LO 

vO 

st 

LO 

o 

r^ 

CD 

t^ 

to 

CD 

CD 

r-- 

CD 

LO 

O 

rH 

aj 

P 

CD 

aj 

r-^ 

tO 

cd 

CX) 

CD 

LO 

St 

CD 

rH 

t-^ 

CN] 

CO 

to 

to 

CN) 

CO 

CD 

CD 

st 

CN 

st 

to 

t-- 

CD 

■p 

g 

r* 

C! 

fH 

OS 

3 

rH 

LO 

CO 

CO 

to 

St 

"vj- 

Tf 

LO 

vO 

r~- 

CO 

CN] 

CN) 

CN) 

rH 

vD 

CNI 

CN 

st 

CD 

st 

o 

vO 

O 

O 

rH 

1 

i— 1 

CO 

t-- 

CD 

LO 

vo 

st 

LO 

CT> 

vO 

r--: 

to 

si- 

o 

o 

C. 

O 

CO 

vO 

rH 

CO 

CO 

CD 

st 

to 

<-H 

crj 

r- 

st 

co 

LO 

vG 

St 

st 

to 

r^ 

rH 

rH 

LO 

to 

vO 

r-i 

CN 

•H 

^ 

vO 
rH 

rH 

rH 
LO 

CO 
rH 

rH 

fH 

o 

Ci 

o 

o 

J> 

A 

•H 

u 

•P 

m 

d> 

aj 

3 

x: 

•H 

cr 

CJ 

■p 

> 

do 

aj 

CN] 

to 

CN 

O 

LO 

CN] 

O 

CD 

to 

c. 

o 

-r 

o- 

O 

st 

O 

to 

CN] 

CD 

VO 

vO 

O 

'ODj 

CN 

CD 

i 

H 

CO 

LO 

r^. 

st 

to 

to 

to 

rH 

r-- 

o 

CD 

VO 

to 

o 

CN) 

i-i 

rH 

CD 

CO 

r-. 

LO 

LO 

st 

LO 

fH 

T3 

H 

01 

o 

_cj 

to 

CD 

00 

CD 

CD 

LO 

CN] 

vO 

LO 

LO 

st 

r-- 

o 

t^ 

CN] 

CO 

LO 

vO 

CO 

to 

c 

CO 

CO 

VO 

vO 

m 

NDJ 

U 

3 

CO 

o 

vO 

O 

to 

CD 

r>- 

CO 

CD 

t^ 

o- 

st 

o 

rH 

to 

o 

O 

LO 

oo 

CN] 

rH 

CO 

CD 

CO 

rH 

aj 

CO 

CNI 

co 

rH 

CO 

st 

CD 

rH 

to 

to 

rH 

st 

to 

to 

CN 

rH 

st 

to 

CD 

aj 

CJ 

rH 

CN] 

rH 

rH 

rH 

u 

13 

CJ 

CJ 

aj 

aj 

•H 

•  -P 

to 

CO 

tO 

o 

I-- 

to 

O 

CD 

r^- 

to 

OD 

co 

O 

LO 

to 

O 

\0 

CD 

CD 

OD 

>-C" 

vO 

rH 

st 

CN 

fn  T^    (A 

•H 

VD 

to 

CO 

CO 

LO 

CO 

st 

st 

CD 

r~) 

o 

rH 

rH 

r^ 

CN) 

st 

rH 

r-i 

vO 

o 

CN 

rH 

to 

CD 

aj 

CD 

W 

CD 

to 

CXI 

CD 

LO 

to 

O 

CD 

\D 

rH 

vO 

CO 

CD 

CD 

to 

to 

CD 

st 

to 

LO 

to 

st 

rH 

CO 

> 

fH     •« 

o 

cd 

rH 

CN 

vO 

rH 

st 

LO 

CD 

CD 

CN] 

st 

Tf 

rH 

to 

VO 

o 

st 

vO 

CO 

to 

LO 

r>- 

CD 

to 

aj   ^) 

•p 

i>~ 

rH 

to 

LO 

LO 

a-, 

F- 

CD 

rH 

O 

CD 

o 

to 

to 

LO 

to 

st 

rH 

st 

vH 

to 

vO 

LO 

r-^ 

.  »> 

2     fH        • 

CO 

u 

t3 

cr  cu  13 

CJ 

3 

CN 

to 

to 

to 

LO 

CN] 

vD 

CN) 

st 

to 

st 

to 

to 

CO 

CN]  K 

st 

st 

to 

to 

to 

LO 

st 

to 

CD 

t/>  -P    C! 

1— 1 

w 

to 

rH 

fH 

CD    O 

E-h 

aJ  tj  S   U 

CO 

J3 

C!          <D 

1— 1 

cr  o  cJ  <s) 

H 

to 

vO 

rH 

OD 

LO 

LO 

CN] 

\o 

to 

VO 

cr, 

t^ 

to 

Tj- 

to 

00 

LO 

CD 

CN) 

oo 

to 

CN 

CO 

vO 

vO 

to 

U  -H 

< 

to 

r-. 

rH 

LO 

vO 

to 

CO 

o 

CO 

o 

CO 

CO 

c  ) 

tO 

CD 

st 

CD 

to 

st 

to 

CN 

st 

t>- 

CO 

to 

CD           fn 

H 

CD 

tO 

t-- 

LO 

to 

st 

LO 

tO 

LO 

CO 

CD 

CD 

CN) 

r-- 

LO 

r>- 

r-~ 

CD 

to 

r>. 

CD 

vO 

CN 

CD 

st 

to 

CD    C    CD 

CO 

£ 

vO 

to 

CD 

h>- 

LO 

I>~ 

o 

CD 

CD 

O 

r-- 

r>- 

rH 

t^- 

co 

CD 

to 

E-» 

r-- 

vO 

t~» 

LO 

to 

st 

fH 

CD    P_ 

£ 

vD 

to 

OD 

(XI 

CD 

rH 

OD 

rH 

st 

rH 

vO 

o 

O 

CD 

o 

LO 

CD 

st 

o 

l>~ 

r-{ 

LO 

CN 

VO 

CD 

fH     > 

n 

CD 

•P 

CD  -H    to 

w 

^! 

o 

o 

O 

CD 

to 

CD 

o 

CD 

rH 

O 

o 

o 

O 

O 

o 

CD' 

CD 

O 

o 

o 

CD 

o 

CD 

O 

CD 

P-,  SO  U 

H 

CO 

1 

I 

1 

I 

1 

e 

CD 

ED 

13    10    -P 

Cm 

CJ 

CD  -H    CD 

O 

CJ 

•H 

fH                  g 

aj   to 

CJ 

13    O 

st 

rH 

LO 

rH 

to 

St 

o 

CD 

CO 

r--. 

CD 

LO 

st 

to 

vO 

CD 

CO 

st 

vO  CO 

vO 

o 

LO 

CN 

L-i 

2  -P    CJ 

fn  -H 

r>- 

rH 

rH 

O 

rH 

rH 

o 

Tj- 

st 

LO 

co.' 

CO 

vO 

CD 

!>- 

H 

CN] 

st 

rH 

CD 

vO 

o 

CO 

rH 

CD 

CTrd     -H 

aj  -P 

CO 

to 

CD 

CM 

CD 

LO 

LO 

tO 

o 

rH 

vO 

t>- 

CD 

r-i 

to 

st 

1^ 

rH 

t>. 

st 

vD 

CD 

to 

VO 

> 

to    bO 

• 

t3   aS 

rH 

VO 

O 

tO 

st 

O 

vO  vO 

t~- 

vO 

CNI 

CN] 

CN) 

CN] 

CN] 

to 

CO 

st 

CN 

i— 1 

CD 

CD 

CD 

CD 

•H 

•h  a 

r-H 

Ci-H 

Osl 

cd 

rH 

CO 

rH 

vO 

rH 

\o 

CD 

i-O 

rH 

st 

rH 

rH 

st 

rH 

LO 

rH 

LO 

rH 

LO 

t-i 

vO 

o 

bfi 

to    CD    CD 

1— 1 

03    > 

fH^> 

+->  a> 

o 

o 

O 

CD 

CD 

o 

CD 

o 

O 

o 

O 

o 

O 

O 

O 

o 

o 

O 

O 

O 

O 

O 

o 

o 

to 

CD          -H 

w 

CO  Q 

•H 

•P    CD    W) 

H-H1 

CD    > 

PQ 

to 

g   aj   to 

<C 

•P 

r5  -H 

H 

cu 

rC 

Ci 

u 

LO 

CO 

LO 

Tf 

st 

o 

to 

CN] 

r-- 

VO 

o 

LO 

rH 

t^- 

CO 

rH 

o 

rH 

CD 

CN 

rH 

O 

to 

st 

bfl 

•H    CD    to 

C! 

CO 

LO 

vo 

CNI 

CN] 

CN] 

CN) 

CN! 

CD 

to 

rH 

CD- 

CO 

CO 

LO 

O 

CD 

CD 

OD 

CN 

O 

vO 

st 

CN 

•H 

r^      +J 

aj 

r-^ 

r>- 

LO 

CN] 

CN] 

VO 

t>- 

CD 

st 

LO 

vO 

to 

vO 

Tf- 

CD 

CO 

st 

O 

r^ 

to 

vO 

CD 

rH 

CD 

CD 

CJ    -P     CD; 

•  H 

st 

CN 

CN 

CD 

CN] 

vO 

CN] 

st 

CD 

rH 

rH 

00 

r-{ 

rH 

r-- 

rH 

st 

CN] 

r-~ 

rH 

LO 

to 

oo 

O 

X 

CD           CD 

fH 

o 

cd 

O 

VD 

CD 

to 

o 

st 

O 

to 

O 

rH 

CD 

O 

rH 

O 

to 

o 

C  ] 

O 

to 

CD 

st 

o 

>  h  rt 

aj 

CD 

•H    O    O 

> 

o 

o 

O 

O 

o 

o 

o 

CD 

CD 

o 

o 

o 

O 

o 

O 

o 

O 

o 

o 

o 

O 

CD 

o 

o 

> 

aj 

IS  g 

on  £ 
comp 

-P 

CD 

•H 

C! 

rC 

to   -P    X 

cd 

-P 

•P   aj  +-> 

6 

rH 

rH 

cn 

CN] 

rH 

rH 

CN] 

CNI 

rH 

rH 

CN) 

CN) 

to 

rH 

rH 

OJ 

CN] 

rH 

rH 

CN 

CN 

rH 

rH 

CN 

C!  -H  -H 

3 

=**= 

=8= 

=* 

3fc 

3fc 

=tfc 

=tb 

=tfe 

4fc 

=*te 

*= 

=«= 

*= 

* 

=tfc 

*: 

3fc 

3fc 

*: 

* 

=te 

=*fc 

=tfe 

=«= 

fH 

CD    >    U 

u 

O 

CJ    CD    O 

-p 

cj 

CJ 

2 

CJ 

2 

CJS 

CJ   S 

S 

CJ 

CJ  S  cJ 

2 

CJ 

CJ  2 

CJ 

P 

CJ 

m 

O  13  rH 

to 

^  P-i 

^ 

Ph 

'^2 

Ph 

^ 

Ph   ^  P-, 

& 

^ 

Pi 

^ 

CL, 

jg: 

Ph 

tg 

Ph 

&   P 

^ 

P_          CD 

cJ 

CD 

S  13   > 

HH 

r] 
U 

o 

u 

rH 

U 

u 

o 

>> 

aj 

u 

C! 

aj 

•H 

fH 

aj 

r> 

city  co 
Standar 
zontal 

CD 

rt 

aj 

crj 

rt 

crj 

*5>* 

O           -H 

■P 

^ 

S 

s 

S 

rH             fH 

aj 

CD 

CD           O 

C 

3 

st 

LO 

VO 

co 

CD 

CN 

>               rC 

28 


s 


CD 

o3 


> 

03 


^ 

v — ' 

v— ' 

O 

h— i 

LO 

LO 

CTi 

o 

rg 

(NI 

H 

VO 

rH 

LO 

CO 

r-» 

en 

!=> 

LO 

O 

LO 

(Nl 

• 

• 

^ 

PQ 

r-g 

vD 

CTl 

eg 

CTl 

CO 

o 

i— i 

o 

i— 1 

CO 

bO 

vO 

t^ 

r— 1 

C* 

• 

• 

• 

• 

CO 

CO 

H 

H 

o 

O 

o 

to 

rH 

D 

00 

PQ 

t— I 

1 

I 

1 

i 

1 

1 

t— i 

Q 

w 

rt 

1 

1 

1 

i 

e*  e* 

I 

1 

H 

& 

<  w 

&, 

to 

m 

1 

I 

1 

i 

D  H 

I 

w 

1— 1 

t— i 

cyw 

r— 1 

Q 

Pi 

! 

I 

1 

i 

oo  s 

I 

►J 

^ 

< 

m 

Q  O 

00 

00 

r-H   K 

Z 

< 

< 

tc 

u 

Oi  t-n 

00 

1— 1 

a:  < 

< 

DC 

y~i 

u 

z 

<  H 

w 

00 

u  a, 

1—1 

u 

00 

i 

< 

Q  < 

2: 

o 

00 

i 

CO 

s 

i— i 

2    t-H 

^ 

E-i 

00 

g 

E> 

2 

C* 

<    > 

w 

P< 

D 

<£ 

< 

k 

< 

H  W 

« 

£3 

< 

K 

LD 

U 

> 

00  Q 

00 

w 

u 

O 

j7 

o 


•p 

o3 

-C 
O 
•H 
■P 

03 
> 
0 


u 

o 

0) 

o 

to 

If) 

^: 

*eT 

£ 

o 

^— ^ 

» — ' 

t—i 

H 

"3- 

i— i 

LO 

o. 

P 

<N) 

o 

vO 

CO 

^ 

PQ 

CN1 

CN) 

tO 

CD 

o 

LO 

o 

►— < 

CT> 

to 

r-. 

^O 

«* 

C-J 

(-H 

Pi 

^O 

CO 

r^ 

LO 

• 

• 

H 

H 

• 

• 

• 

• 

cr> 

CO 

n 

CO 

O 

o 

o 

to 

O 

Ol 

PQ 

h- 1 

w 

rH 

h— i 

Q 

I 

1 

i 

i 

p:  pi 

I 

I 

Pi 

<  w 

H 

Pi 

I 

1 

i 

i 

ZD  H 

I 

Pi 

CO 

W 

C/W 

w 

t— i 

t— i 

I 

1 

i 

i 

CO  s 

I 

K- 1 

O 

Pi 

, 

^ 

i 

i 

>— 1  p. 

, 

—1 
Pi 

2 

< 

w 

Q  O 

CO 

CO 

a:  < 

•z 

< 

< 

E 

u 

Pi   1-H 

CO 

1— < 

U  P-, 

< 

X 

i— i 

CJ 

:s 

<  H 

w 

CO 

1— 1 

u 

CO 

1 

< 

Q< 

z 

o 

CO 

1 

CO 

X 

H-l 

2:  m 

V: 

H 

CO 

T 

;=» 

^ 

<  > 

W 

Pi 

D 

< 

< 

c2 

< 

H  W 

X 

ED 

< 

Pi 

o 

C_3 

> 

CO  Q 

CO 

^ 

o 

u 

rsj 

4t= 

P 

O 

+-> 

CD 

S 

£ 

o 

i-H 

Ph 

+-> 

03 

>n 

•P 

•H 

O 

O 

i— 1 

CD 

> 

i-H 

nJ 

+->    LO 

£    ^ 

O  cr> 

N  r- 1 

•  H 

mx: 

o  o 

K    ^ 

0J 

<L»   S 

X 

■P    <^ 

CO 

2 

m  c 

o 

o  o 

t— 1 

H 

c 

< 

o 

i—i 

•H 

> 

•P 

W 

3 

Q 

rO 

•H 

O 

h 

Pi 

p 

< 

w 

Q 

•H 

^ 

Q 

< 

H 

>s 

CO 

U 

£ 

0 

3 

cr 

CD 

P 

Ph 

w 
Pi 
:=> 

i— i 

Ph 


30 


wave  components.   This  is  to  be  expected  since  the  flow 
meters  do  not  experience  the  surface  irregularities  to 
the  same  degree  as  do  the  wave  gauges.   In  seven  of 
thirteen  fit  parameter  calculations  for  wave  components, 
the  Gaussian  distribution  provided  the  better  fit.   This 
is  a  result  of  using  a  truncated  form  of  the  Gram-Charlier 
in  calculating  the  probability  density  function,  Pfr , 
given  by 

m,      m.      m_ 

pGCco  -  pGcmi.o  ♦  ^Hj  ♦2|h4  ♦  jf-Hg  ♦  ]    c 


where 

P„  =  Gaussian  probability  density  function, 

m~  =  skewness  of  the  data  record, 

m.  =  the  kurtosis  minus  3, 

H   =  Hermite  polynomials  of  degree  n. 
Higher  order  moments  could  be  included  in  equation  (1) 
to  possibly  improve  the  fit,  but  would  necessitate  analyzing 
longer  record  lengths  in  order  to  maintain  confidence. 
Stationary  problems  would  be  encountered  in  the  longer 
records,  thus  producing  other  errors.   Hence,  the  poor 
fit  is  a  deficiency  of  the  Gram-Charlier  distribution  in 
surf  zone  applications. 

2 .   Skewness  and  Kurtosis 

The  physical  geometry  of  the  waves,  which 
characteristically  have  narrow,  steep  crests  and  wide, 
shallow  troughs  in  the  surf  zone  produce  a  positive 
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skewness,  indicating  a  greater  amount  of  time  below, 
rather  than  above,  the  mean  water  level.   This  is  true 
except  at  Wave  Gauge  #3  on  6  March.   (See  Table  II  and 
Appendix  F.)   Data  from  this  wave  gauge  was  analyzed  on 
this  date  only,  because  visual  observations  noted  that 
most  waves  were  actually  breaking  on  this  gauge,  rather 
than  on  Wave  Gauges  #1  and  #2  as  was  otherwise  the  case. 

It  should  follow  that  the  wave  induced  particle 
velocities  would  have  probability  density  functions  similar 
to  the  waves.   However,  of  the  eleven  horizontal  velocity 
data  records,  four  did  not  have  positive  skewness.   These 
deviations  occurred  at  Flow  Meter  #1  on  4  March  and 
29  May  and  Flow  Meter  #2  on  5  and  6  March.   Reflected 
waves  were  originally  thought  to  be  a  feasible  explanation 
for  the  negative  skewness.   If  this  were  so,  Flow  Meter  #2 
on  A    March  and  29  May  would  also  expect  to  have  negative 
values.   This  was  not  the  case.   The  reason  for  the  large 
skewness  and  kurtosis  values  for  Wave  Gauge  #1  on  5,  6  and 
8  March  is  not  known.   Similar  anomalous  results  have  been 
found  in  deeper  water  outside  the  surf  zone  by  Thornton 
and  Krapohl  (1974) . 

Kurtosis  indicates  the  peakedness  of  each 
parameter.   Just  prior  to  breaking  the  waves  achieve  the 
greatest  degree  of  peakedness  and  should  have  the  highest 
kurtosis  value  at  this  time.   Visual  observations  noted 
that  most  waves  were  breaking  at  Wave  Gauge  HI    on 
4  March  and  29  May,  at  Wave  Gauge  #2  on  5  March  and  on 
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Wave  Gauge  #3  on  6  March.   It  should  be  expected  that  the 
kurtosis  values  would  be  greatest  at  these  given  gauges 
for  the  respective  dates.   With  the  exception  of  29  May, 
this  did  not  turn  out  to  be  true.   The  wide  surf  zone 
caused  by  varying  breakers  heights  on  8  and  9  March  make 
evaluation  of  kurtosis  values  difficult  on  these  dates. 

C.   COMPARISONS  OF  THEORETICAL  AND  MEASURED  POWER  SPECTRA 

The  theoretical  velocity  spectra  generated  from  the 
wave  spectra  and  the  measured  velocity  spectra  were 
examined  in  order  to  determine  the  applicability  of  linear 
(Airy)  wave  theory  in  the  breaker  zone.   Thornton  and 
Richardson  (1974)  measured  waves  with  maximum  heights  of 
1  m  and  showed  that  the  application  of  a  linear  theory 
transfer  function  to  the  surface  profile  spectrum  resulted 
in  a  calculation  of  the  horizontal  water  particle  velocity 
that  underpredicted  measured  wave-induced  horizontal 
velocity  spectral  components  by  about  50  percent.   The 
coherence  values  between  waves  and  horizontal  velocity 
were  high,  ranging  above  0.75  for  an  approximate  frequency 
range  of  0.075  to  0.60  Hz.   The  phase  angle  computations 
showed  the  calculated  velocity  components  leading  the 
measured  velocity  components  by  an  average  of  20  degrees. 
Bub  (1974)  took  measurements  in  mild  surf  resulting  in 
theoretical  horizontal  velocity  values  13  percent  lower 
than  measured  velocity.   Coherence  values  varied  from 
0.5  to  nearly  1.0  over  a  frequency  range  of  0.10  to  0.65  Hz 
An  approximate  zero  degree  phase  was  observed  out  to  1.0  Hz 

33 


For  this  study  wave  heights  of  approximately  1  m 
were  typically  observed  at  Del  Monte  Beach  and  up  to  2  m 
were  recorded  at  Carmel  River  Beach.   The  theoretical 
horizontal  velocity  spectra  were  computed  from  the  wave 
spectra  using  the  linear  theory  transfer  function 

S   m  -  \°   cosh(k(h+z))-.2    rr,  r?. 

ucl±J  -  l sinh(kh) ]   Vf)  (2) 


where 


S   (f)  =  calculated  horizontal  velocity  power 
uc  }    l 

2 
spectrum  (m  /sec)  , 
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S  (f)   =  measured  wave  height  power  spectrum  (m  -sec), 

a  =  angular  frequency  =  2TTf  (sec   )  , 

k  =  wave  number  (m   ) , 

h  =  mean  water  depth  (m) , 

z  =  depth  of  flow  meter  below  the  mean  water 
depth  (m) . 
Figures  8  and  9  are  graphs  of  the  measured  and  calculated 
spectral  velocity  values  at  Tower  #1  for  4  March  and 
29  May.   Additional  spectra  are  contained  in  Appendix  G. 
On  the  average  at  Del  Monte  Beach,  the  calculated  hori- 
zontal velocity  values  underestimated  the  measured  veloci- 
ties by  7  9%.   The  same  theoretical  calculations  for  the 
data  collected  at  Carmel  River  Beach  underestimated  the 
measured  values  by  86%. 

A  possible  explanation  for  the  low  calculated  values 
is  that  the  wave  profiles  change  during  the  shoaling  and 
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breaking  process,  becoming  more  peaked  at  the  crest, 
thus  lessening  the  applicability  of  linear  (Airy)  theory 
which  assumes  sinusoidal  wave  shapes.   Additionally,  the 
wider  separation  between  the  measured  and  theoretical 
spectral  values  above  0 . 5  Hz  may  be  explained  by  turbu- 
lence of  the  breaking  waves  in  the  higher  frequency  areas . 

D.   EXAMINATION  OF  SPECTRAL  HARMONICS 

An  examination  of  the  low  frequency  spectral  components 
for  each  data  set  indicates  that  pronounced  approximate 
harmonics  are  present  in  the  data  analyzed  from  Del  Monte 
Beach.   Such  harmonics  are  not  evident  in  the  Carmel  River 
Beach  spectra.   The  frequencies  and  amplitudes  of  each 
analyzed  spectral  peak  is  listed  in  Table  III.   The  fre- 
quency differences  between  spectral  peaks  and  the  ratios 
of  the  second,  third  and  fourth  peaks  versus  the  primary 
peak  are  listed  in  Table  IV. 

Spectral  analysis  shows  that  definite  second,  third 
and  fourth  harmonics  are  present  for  4  March  and  5  March 
with  few  exceptions .   Evidence  of  second  and  third  har- 
monics is  present  in  the  analyzed  data  from  Wage  Gauges 
#1  and  #2  on  6,  8  and  9  March.   An  examination  of  the 
spectral  amplitudes  shows  that  there  is  no  consistent 
relationship  between  various  components  either  for  a 
given  data  set  or  when  comparing  like  components  of 
different  data  sets. 

The  maximum  wave  energies  were  concentrated  at  frequen- 
cies ranging  from  0.061  Hz  to  0.092  Hz  corresponding  to 
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periods  of  16.4  to  10.9  seconds,  respectively,  which  agreed 
with  visual  observations.   Additionally,  prominent  sub- 
harmonic  peaks  between  0.011  and  0.022  Hz  are  noted  in  the 
measured  velocity  spectra  computed  from  Del  Monte  Beach 
data.   This  agrees  with  the  work  of  Bub  (1974),  who 
observed  a  sub-harmonic  at  0.011  Hz.   Guzza  and  Davis  (1974) 
state  that  in  a  theoretical  analysis  of  edge  waves  excited 
by  incoming  waves,  the  prominent  edge  wave  mode  is  the 
first  sub-harmonic  of  the  primary  frequency.   It  is  plausi- 
ble, therefore,  that  the  observed  low  frequency  sub -harmonic 
peaks  were  caused  by  edge  waves  in  the  surf  zone.   Such 
sub-harmonics  are  not  observed  in  the  measured  velocity 
spectra  generated  from  Carmel  River  Beach  data. 

E.   PHASE  AND  COHERENCE 

The  curling  crests  of  the  unstable  breaking  waves 
tend  to  lead  the  horizontal  particle  velocities  in  the  body 
of  the  waves  as  can  be  seen  in  the  phase  measurements 
between  waves  and  velocities.   (See  Table  V.)   The  maximum 
phase  difference  would  be  expected  at  breaking.   Data  from 
Tower  #1  on  4  March  indicates  that  the  theoretical  hori- 
zontal velocity  leads  the  measured  horizontal  velocity  by 
about  20  degrees.   Since  most  waves  observed  were  breaking 
on  Tower  #1,  a  possible  explanation  for  the  30  degree  phase 
angle  at  Tower  2  is  that  the  waves  reorganized  between  the 
towers  and  broke  a  second  time  at  Tower  #2.   Since  this 
author's  observations  concentrated  primarily  on  the  initial 
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breaking  position,  the  possible  reorganization  and  second 
breaking  were  not  observed.   The  phase  angles  on  5  March 
are  5  and  20  degrees  at  Towers  #1  and  #2,  respectively. 
This  agrees  with  observations  and  indicates  that  the 
breakers  started  to  curl  at  Tower  #1  and  broke  at  Tower  #2. 
The  zero  degree  phase  relationships  on  6  March  show  that 
the  waves  had  not  started  to  break  until  shoreward  of 
Tower  #2.   This  agrees  with  the  observation  of  waves 
breaking  on  Tower  #3.   The  phase  angles  on  8  and  9  March 
indicate  that  the  smaller  breakers  began  to  curl  at  Tower  #1 
and  were  near  breaking  at  Tower  #2.   Again,  this  agrees 
with  observations. 

In  contrast  to  the  plunging- type  breakers  at  Del  Monte 
Beach,  surging- spilling  breakers  were  observed  at  Carmel 
River  Beach.   This  was  to  be  expected  considering  the  beach 
slopes.   A  spilling  breaker  is  characterized  by  turbulent 
water  forming  at  the  wave  crest  and  eventually  flowing  down 
and  covering  the  leading  edge  of  the  wave.   Spilling  begins 
at  the  crest  when  a  small  tongue  of  water  moves  faster 
than  the  wave  form  as  a  whole  (Galvin,  1972).   The  measured 
phase  angle  at  Carmel  River  Beach  shows  that  the  cascading 
crests  led  the  horizontal  particle  velocities  in  the  body 
of  the  wave  by  10  degrees. 

The  values  of  coherence  show  varying  degrees  of  linearity 
in  the  transfer  functions  between  different  spectra.   Co- 
herence between  wave  and  horizontal  velocity  spectra  at 
Tower  #1  on  29  May  is  above  0.53  out  to  0.S8  Hz  (See 


4  2 


Figure  9),  above  0.68  out  to  0.53  Hz  on  8  March  and  above 
0.60  out  to  0.41  Hz  on  9  March.   (See  Appendix  G.)   Other 
coherence  calculations  show  a  fair  to  poor  degree  of 
linearity,  which  is  consistent  with  the  observed  strong 
non- linearities . 

F.   THEORETICAL  AND  CALCULATED  PHASE  SPECTRA 

A  theoretical  phase  spectrum  was  calculated  for  com- 
parison with  the  measured  wave  spectrum  from  two  wave 
gauges.   The  theoretical  phase  was  computed  from  the 
relationship 

*  =  P  '         (3) 

Ll 

where 

4>  =  phase, 

a  =  angular  frequency  =  2irf  (sec   )  , 

x  =  distance  between  wave  gauges  (m) , 

1  /? 
C,  =  celerity  =  (gh)  '   (m/sec) , 

2 
g  =  gravity  (m/sec  )  , 

h  =  mean  water  depth  (m) . 

Figure  10  shows  typical  examples  of  the  theoretical  and 

measured  phases.   The  theoretical  phases  computed  from  Del 

Monte  Beach  data  underpredict  the  actual  phase  relationships 

whereas  an  overprediction  is  noted  in  the  Carmel  River  Beach 

phases.   The  author  then  determined  what  value  of  celerity 

(C?)  would  cause  the  theoretical  and  actual  phases  to 

agree  out  to  the  first  zero  crossing  (where  <j)  =  180°  =  tt) 
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FIGURE  10.   Examples  of  Theoretical  and  Measured  Phases 
from  Del  Monte  Beach  (top)  and  Carmel  River 

Beach 
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by  first  solving  for  C?  from  an  iteration  of  equation  (3) 


C2  =  2fx  (4) 


and  then  solving 


C2  =  (g(h+a))1/2  (5) 

for  a.   The  results  are  shown  in  Table  VI.   No  evident 
relationship  between  a  and  other  parameters,  such  as  beach 
slope,  mean  water  depth  and  standard  deviation,  could  be 
determined. 
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V.   CONCLUSIONS 

Breaking  waves  in  the  surf  zone  can  be  characterized 
as  being  highly  non-linear.   Qualitative  observations  of 
wave  and  velocity  profiles  show  secondary  waves  indicative 
of  non-linear  waves.   The  probability  density  functions 
calculated  for  the  flow  velocity  records  compared  better 
with  the  Gram-Charlier  distribution  than  with  the  Gaussian 
distribution  when  tested  using  the  chi-square  goodness-of- 
fit  test.   However,  seven  of  the  thirteen  wave  records 
more  closely  approximated  the  Gaussian  distribution.   This 
is  a  result  of  the  truncated  form  of  the  Gram-Charlier 
probability  density  function  used  in  calculations,  and 
points  out  the  importance  of  including  higher  order  moments 
in  describing  wave  phenomena  in  very  shallow  water. 

The  values  of  the  horizontal  power  spectral  components 
calculated  from  wave  spectra  using  linear  theory  indicate 
a  qualitative,  but  not  a  quantitative,  relationship.   Com- 
paring the  results  of  Bub  (1974) ,  Thornton  and  Richardson 
(1974)  and  this  author,  it  is  evident  that  increased 
breaker  heights  result  in  larger  estimation  errors. 
Additional  investigation  is  warranted  to  determine  if 
specific  breaker  heights  and/or  types  result  in  specific 
underestimation  values.   If  so,  linear  theory  may  be 
applied  to  achieve  a  reasonably  accurate  determination  of 
flow  velocities.   The  coherence  values  between  waves  and 
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horizontal  velocity  decrease  as  breaker  height  increases, 
indicating  that  the  wave  motion  becomes  more  non-linear 
and  more  turbulent  under  this  condition. 

The  theoretical  phase  calculations  based  on  the  linear 
theory  approximation  for  celerity  did  not  accurately 
predict  the  observed  phases  in  this  research.   The  values 
of  a,  determined  so  as  to  force  the  theoretical  phases  to 
agree  with  the  observed  values  at  the  first  zero  crossing, 
show  no  consistency. 

Use  of  electromagnetic  flow  meters  and  the  improved 
capacitance  wave  gauges  and  instrument  towers  permits  the 
gathering  of  accurate  continuous  data  with  sturdy  and 
reliable  equipment. 
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APPENDIX  A 
BEACH  PROFILES  AT  DEL  MONTE  BEACH  4-10  MARCH  1975 
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APPENDIX  B 
CAPACITANCE  WAVE  GAUGE  CALIBRATIONS 


x  r  wave  gauge  1 
a  z  wave  gauge  2 
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APPENDIX  C 
FLOW  METER  CALIBRATIONS 
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APPENDIX  D 
ANALYSIS  DETAILS  AND  FLOW  CHART 

A  complete  flow  chart  of  data  acquisition  and  analysis 
is  presented  on  the  following  page.   The  parameters  used 
during  spectral  analysis  were  chosen  to  balance  the  longest 
record  that  the  computer  could  reasonably  analyze  with  the 
best  resolution  over  the  frequency  range  of  interest. 
Additionally,  consideration  was  given  to  the  computer  run 
time . 

In  order  to  determine  the  optimum  Nyquist  frequency  for 
data  analysis  the  following  equations  were  used: 

At  =  (VDT) (NSKIP) (NCHAN)  (1) 

and 

fN  =  TfJtJ  (2) 

where 

At  =  sampling  interval, 
VDT  =  sampling  rate, 
NSKIP  =  number  of  samples  skipped  in 

initial  data  array, 
NCHAN  =  number  of  channels, 
f   =  Nyquist  frequency. 
By  varying  the  parameters  of  equation  (1)  various  Nyquist 
frequencies  were  obtained  and  used  during  preliminary  IBM 
360  computer  analysis. 
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APPENDIX  E 
CALIBRATION  FACTORS 


GALA*  CALM** 

(meters)        (meters/volt) 


Date  Instrument 


4  March 

WG  #1 

and 

WG  #2 

5  March 

WG  #3 

6  March 

WG  #1 

WG  #2 

WG  #3 

8  March 

WG  #1 

WG  #2 

WG  #3 

9  March 

WG  #1 

WG  #2 

WG  #3 

29  May 

WG  #1 

WG  #2 

All  dates 

FM  #1 

FM  #2 

-1.312  1.23 

-1.870  4.19 

-1.800  1.14 

-1.422  1.23 

-1.620  4.19 

-1.800  1.14 

-1.332  1.23 

-1.570  4.19 

-1.810  1.14 

-1.492  1.23 

-1.510  4.19 

-1.830  1.14 

-0.028  0.69 

-0.370  0.66 

0.0  3.05 

0.0  3.21 


*  Calibration  Additive  Factor 
**  Calibration  Multiplicative  Factor 
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APPENDIX  F 

PROBABILITY  DENSITY  FUNCTIONS  VS. 
GRAM-CHARLIER  FREQUENCY  DISTRIBUTIONS 
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APPENDIX  G 
POWER,  COHERENCE  AND  PHASE  SPECTRA 

FLOW  METER  #2  -  WAVE  GAUGE  #2 
4  MARCH  1975 
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FLOW  METER  #1  -  WAVE  GAUGE  #1 
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